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Abstract 

Merging junctions are important network bottlenecks, and a better understanding of merg- 
ing traffic dynamics has both theoretical and practical implications. In this paper, we present 
continuous kinematic wave models of merging traffic flow which are consistent with discrete 
Cell Transmission Models with various distribution schemes. In particular, we develop a sys- 
tematic approach to constructing kinematic wave solutions to the Riemann problem of merging 
traffic flow in supply-demand space. In the new framework, Riemann solutions on a link con- 
sist of an interior state and a stationary state, subject to admissible conditions such that there 
are no positive and negative kinematic waves on the upstream and downstream links respec- 
tively. In addition, various distribution schemes in Cell Transmission Models are considered 
entropy conditions. In the proposed analytical framework, we prove that the stationary states 
and boundary fluxes exist and are unique for the Riemann problem for both fair and constant 
distribution schemes. We also discuss two types of invariant merge models, in which local 
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and discrete boundary fluxes are the same as global and continuous ones. With numerical 
examples, we demonstrate the validity of the analytical solutions of interior states, stationary 
states, and corresponding kinematic waves. Discussions and future studies are presented in the 
conclusion section. 

Key words: Kinematic wave models; merging traffic; Riemann problems; Cell Transmission Mod- 
els; distribution schemes; supply-demand space; stationary states; interior states; boundary fluxes; 
invariant merge models 



1 Introduction 



In a road network, bottlenecks around merging, diverging, and other network junctions can cause 
the formation, propagation, and dissipation of traffic queues. For example, upstream queues can 
form at a merging junction due to the limited capacity of the downstream branch; at a diverg- 
ing junction, when a queue forms on one downstream branch, the flow to the other downstream 
branch will be reduced, sin ce the upstream link will be blocked due to the First-In-First-Out prin- 



ciple on the upstream link ( 



Papageorgioul . 



1990h . Theoretically, a better understanding of merging 



traffic flow wil 



( Nie andZhand. 



DC helpfu 



2008 



Jin, 



for an alyzing and understanding traffic dynamics on a road network 



(|Papageorgiou and Kotsialos 



20091). Practically, it will be helpful for designing metering schemes 



20021) and for understandi ng dynamic user equilibrium and drivers' 



route choice behaviors (|Peeta and Ziliaskopoulos , 



20011) 



In the liter ature, a number of m icroscopic models have been proposed for vehicles' merging 



behaviors (e.g. 



HidasL 



2001 



2005 



. Due to the complicated interactions among vehicles, micro- 



scopic models are not suitable for analyzing merging traffic dynamics at the network level. At 



the macroscopic 



evel, many attempts have been made to model merging traffic flow in the line of 



the LWR model (ILighthill and Whitham . 



1955 



Richards 



1956h . which describes traffic dynamics 



with kinematic waves, including shock and rarefaction waves, in density (p), speed (v), and flux 



iq). Based on a continuous version of traffic conservation, ^ + ^ = 0> and an assumption of a 
speed-density relationship, v = y(p), the LWR model can be written as 



|p + |py(p)^o, 



(1) 



which is for a homogeneous road link with time and location independent traffic characteristics, 
such as free flow speed, jam density, capacity, and so on. In general, V{p) is a non-increasing 
function, and Vf = V{0) is the free flow speed. In addition, q = Q{p) = pV{p) is unimodal with 
capacity C = Q{pc), where pc is the critical density. If we denote the jam density by pj, then 
pe [0,pj]. 



In one line. 



Daganzd (Il995|) and 



Lebacqud (119961) extended the Godunov discrete form of the 



LWR model and developed a new framework, referred to as Cell Transmission Models (CTM) 
hereafter, for numerically computing traffic flows through merging, diverging, and general junc- 
tions. In this framework, local traffic demand (i.e., sending flow) and supply (i.e., receiving 
flow) functions are introduced, and boundary fluxes through various types of junctions can be 
written as functions of upstream demands and downstream supplies. In the CTM framework, 
various distr ibution schemes can be employed to uniquely determine merging flows from all up- 



stream links (IDaganzo , 



1995 



Lebacqud . 



19961 : 



Jin and Zhang , 



2003 



Ni and Leonard . 



2005 ). These 



discrete merge models are physically intuitive, and stationary distribu tion patterns of rt i erging 



flows have been observe d and calibrated for various merging junctions (ICassidy and Ahn , 



Bar-Gera and Ahn 



2005 



20091). Such merge models are discrete in nature and can be used to simulate 



macroscopic traffic dynamics efficiently. However, there have been no systematic approaches to 
analyzing kinematic waves arising at a merging junc tion with these mod els. 



In the other line. 



Coclite et al 



(|2005|) attempted to solve a Rie- 



Holden and Risebrd (| 19951) and|< 
mann problem of a highway intersection with m upstream links and n downstream links. In both 
of these studies, all links are homogeneous and have the sar ne speed-density relations, a nd traffic 



dynamics on each link are described by the LWR model. In (|Holden and Risebro , 



19951), the Rie- 
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mann problem is solved b y introducing an entr opy condition that maximizes an objective function 



in all boundary fluxes. In (|Coclite et al. . 



2005|) . the Riemann problem is solved to maximize total 
flux with turning proportions. Both studies were able to describe kinematic waves arising from a 
network intersection but also subject to significant shortc omings: (i) All links are ass umed to have 



the same fundamental diagram in both studies; (ii) In (IHolden and Risebro 



19951). the ent ropy 



conditions are pragmatic and lack of physical interpretations; and (iii) In (fCoclite et al. . 



20051), the 



Riemann problem can only be uniquely solved for junctions with no fewer downstream links, and 
hence their results do not apply to a merging junction. 

In this paper, we are interested in studying continuous kinematic wave models of merging 
traffic flow which are consistent with discrete CTM with various distribution schemes. Here we 
consider a merge network with m > 2 upstream links and one downstream link, as shown in Figure 
[T] In this network, there are m + 1 links, m origin-destination pairs, and m paths. In the continuous 
kinematic wave models of this network, traffic dynamics on each link are described by the LWR 
model (dl), and, in addition, an entrop y condition bas ed on various distribution schemes is used 



to pick out unique physical solutions (lAnsorgd . 



19901) . Note that traffic dynamics on the whole 



network cannot be modeled by either one-dimensional or two-dimensional conservation laws in 
closed forms, since vehicles of different commodities interact with each other on the downstream 
link. Similar to the LWR model, it is not possible to obtain analytical solutions for general initial 
and boundary conditions, and we usually resort to analyzing continuous kinematic wave solutions 
to Riemann problems. Studies on the Riemann problem are helpful for understanding fundamental 
characteristics of the corresponding merge model and constructing numerical solutions. In the 
Riemann problem, all links are homogeneous and infinitely long; for link a = l,---,m+l, we 
assume that its flow-density relation is qa = Qa{pa), critical density pc.a, and its capacity Q; and 
upstream link z = 1 , • ■ ■ , m and downstream link m+1 have constant initial conditions: 



p/, X/ G (-00,0), z=l, 



(2) 



4 



Pm+l(-«m+l,0) = Pm+1, -^ra+l e (0,+oo). 



(3) 



To pick out physical kinematic wave solutions to the Riemann problem, here we use various distri- 
bution schemes of C TM as our entropy condition s, which are different from the optimization-based 



entropy condition in (IHolden and Risebro. . il995i) . In this sense, the resulted models can be consid- 



ered as continuous versions of CTM with various distribution schemes. 

1 ^ 




•m + l 



x„ = 



Figure 1: An illustration of a merge network 



In our study, the upstream links can be mainline freeways or on-ramps with different character- 
istics, and our solutions are physically meaningful and consi stent with discre te CTM of merging 



traffic. In addition, we follow a new framework developed in (|Jin et al. 



2009|) for solving the Rie- 



mann problem for inhomogeneous LWR model at a linear junction. In this framework, we first 
construct Riemann solutions of a merge model in the supply-demand space: on each branch, a 
stationary state will prevail after a long time, and an interior state may exist but does no t take any 



space in the continuous solution and only shows up in one cell in the numerical solutions (|van Leer . 



1984h . After deriving admissible solutions for upstream and downstream stationary and interior 
states, we introduce an entropy condition consistent with various discrete merge models. We then 
prove that stationary states and boundary fluxes are unique for given upstream demands and down- 
stream supply, and interior states may not be unique but inconsequential. Then, the kinematic wave 
on a link is uniquely determined by the corresponding LWR model with the stationary state and 



the initial state as initial states. 



In this study, we do not attempt to develop new simulation models. Rather, we attempt to 
understand kinematic wave solutions of continuous merge models consistent with various CTM. 
The theoretical approach is not meant to replace numerical methods in CTM but to provide insights 
on the formation and dissipation of shock and rarefaction waves for the Riemann problem. In this 
study, we choose to solve the Riemann problem in supply-demand space due to its mathematical 
tra ctability. In fact , the benefits o f work ing in the supply-demand space have been demonstrated 



by 



Daganzo 



19951) and 



Lebacqud ( 



19961) when they developed discrete supply-demand models of 
network vehicular traffic, i.e., CTM. In their studies, the introduction of the concepts of demand 
and supply enabled a simple and intuitive formulation of merge and diverge models. But different 
from these studies, which focused on discrete simulation models, this study further introduces 
supply-demand diagram and discusses analytical solutions of continuous merge models. 

The rest of the paper is organized as follows. In Section 2, we introduce an analytical frame- 
work for solving the kinematic waves of the Riemann problem with jump initial conditions in 
supply-demand space. In particular, we derive traffic conservation conditions, admissible con- 
ditions of stationary and initial states, and additional entropy conditions consistent with discrete 
merge models. In Section 3, we solve stationary states and boundary fluxes for both fair and 
constant merge models. In Section 4, we introduce two invariant merge models. In Section 5, 
we demonstrate the validity of the proposed analytical framework with numerical examples. In 
Section 6, we summarize our findings and present some discussions. 



2 An analytical framework 



For link a = 1 , • • • , m -|- 1, we define the following demand and supp l y func tions with all subscript 



1 I 1 r 



a suppressed (.Engquist and Qsheruil980i : iDaganzo . 



19951: 



Lebacquel 



1996 ) 



D{p) = e(min{p,p,}) 



e(p), if p<pc 



if P > Pc 



= r x{s)Q\s)ds= rmax{Q\s),0}ds (4) 
Jo Jo 

Q{p), if p>pc 



5(p) = e(max{p,pj) 



C, if p < pc 
rP 



= C+ [ {l-xis))Q'is)ds = C+ [ mm{Q'{s),0}ds, (5) 
Jo Jo 

where x{p) equals 1 iff Q'{p) >0 and equals otherwise. 

Unlike many existing studies, in which traffic states are considered in p-q space, we represent 
a traffic state in supply-demand space as U = {D,S). For the demand and supply functions in Q 
and ([5]), we can see that D is non-decreasing with p and S non-increasing. Thus D < C, S < C, 
max{D,5} = C, and flow-rate q{U) = min{D,5}. In addition, D = 5 = C iff traffic is critical; 
D < 5 = C iff traffic is strictly under-critical (SUC); 5 < D = C iff traffic is strictly over-critical 
(SOC). Therefore, state U = {D,S) is under-critical (UC), iff S = C, or equivalently D < S; State 
U = (D, S) is over-critical (OC), iff D = C, or equivalently S < D. 

In Figure [2tb), we draw a supply-demand diagram for the two fundamental diagrams in Figure 
[2ta). On the dashed branch of the supply-demand diagram, traffic is UC and U = {D,C) with 
D < C; on the solid branch, traffic is OC and U = (C, S) with S < C. Compared with the funda- 
mental diagram of a road section, the supply-demand diagram only considers its capacity C and 
criticality, but not other detailed characteristics such as critical density, jam density, or the shape 
of the fundamental diagram. That is, different fundamental diagrams can have the same demand- 
supply diagram, as long as they have the same capacity and are unimodal. However, given a 
demand-supply diagram and its corresponding fundamental diagram, the points are one-to-one 
mapped. 

n 

In supply-demand space, initial conditions in Q and ^ are equivalent to\j 

Ui{xuQ) = iD,,Si), x,G(-oo,0), (6) 

U,n+l{Xm+uO) = (D„,+ i,5,„+l), Xm+l e {0,+oo). (7) 



In this section, / = 1 , • • • , w if not otherwise mentioned. 




Figure 2: Fundamental diagrams and their corresponding supply-demand diagrams 

In the solutions of the Riemann problem with initial conditions (I6l|7]), a shock wave or a rarefac- 
tion wave could initiate on a link from the merging junction x = 0, and traffic states on all links 
become stationary after a long time. That is, in solutions to the Riemann problem, statio nary states 



preva i 



19841 : 



all links afte r a lon g time. At the boundary, there can also exist interior states (|van Leei . 



Bultelle et al 



1998 ). which take infinitesimal spaced. We denote the stationary states on 
upstream link i and downstream link m -|- 1 by Uf' and U^^^, respectively. We denote the interior 
states on links i and m+ 1 by Ui{0^ j) and Um+iiO'^ j), respectively. The structure of Riemann 
solutions on upstream and downstream links are shown in Figure |3l where arrows illustrate the 
directions of possible kinematic waves. Then the kinematic wave on upstream link / is the solution 
of the corresponding LWR model with initial left and right conditions of Ui and i7-^, respectively. 
Similarly, the kinematic wave on downstream link m -|- 1 is the solution of the corresponding LWR 
model with initial left and right conditions of U^^-^ and Um+u respectively. 

We denote <7,_>m+i as the flux from link / to link m + 1 for t > 0. The fluxes are determined 
by the stationary states: the out-flux of link / is qt = q{Uj^), and the in-flux of link m + 1 is 
^In numerical solutions, the interior states exist in one cell. 
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r/,(o-.f) 



C,„+i(0+,i) 
J 



U+ 



(a) 



(b) 



Figure 3: The structure of Riemann solutions: (a) Upstream link /; (b) Downstream link m + 1 



<?m+i = <?(f^„+i)- Furthermore, from traffic conservation at a merging junction, we have 

m 

qi^m+l = qi = q{U^), qm+l = q{U^+i) = l^^(f/r)- 



(8) 



i=l 



2.1 Admissible stationary and interior states 



As observed in (Holden and Risebrd . ll995l : ICoclite et al.Ll2005h . the speed of a kinematic wave on 
an upstream link cannot be positive, and that on a downstream link cannot be negative. We have 
the following admissible conditions on stationary states. 

Theorem 2.1 (Admissible stationary states) For initial conditions in © and ([7|), stationary states 
are admissible if and only if 



{Di,Q) or (Q,5r), 



(9) 



where 5- < , and 



\-,Sm+\) or {D^^^^,Cm+\, 



(10) 



where < S^+i ■ 

The proof is quite straightforward and omitted here. The regions of admissible upstream stationary 
states in both supply-demand and fundamental diagrams are shown in Figure H and the regions of 
admissible downstream stationary states are shown in Figure [51 In Figure |4ta) and (b), the initial 
upstream condition is SUC with D,- < Q = 5,-, and feasible stationary states are (D,-,C/), when 
there is no wave, or {Ci.Sj) with SJ^ < D,-, when a back-traveling shock wave emerges on the 



upstream link /; In Figure Stc) and (d), the initial upstream condition is OC with 5, < Di = Ci, 
and any OC stationary state is feasible, when a back-traveling shock or rarefaction wave emerges 
on the upstream link i. In contrast, in Figure [2a) and (b), the initial downstream condition is 
UC with D,n+\ < S,n+i = Cm^i, and any UC stationary state is feasible, when a forward-traveling 
shock or rarefaction wave emerges on the downstream link m+ 1; In Figure [3c) and (d), the 
initial downstream condition is SOC with 5^+1 < £>m+i = Cm+\, and feasible stationary states 
are (C„,+i,5m+i), when there is no wave, or (D+^j,Cm+i) with D'^^-^ < S,n+i, when a forward- 
traveling shock wave emerges on the downstream link m+l. Here the types of kinematic waves 
and the signs of the wave speeds can be determined in the supply-demand diagram, but the absolute 
values of the wave speeds have to be determined in the fundamental diagram. 



s, 



Ui = u- 




-® 







Pi 



(a) 



(b) 



s, 








(c) 



(d) 



Figure 4: Admissible stationary states for upstream link /: marked by black dots 
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(c) 




(b) 





/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 




I 

• 





Pi 



(d) 



Figure 5: Admissible stationary states for downstream link m+ 1: marked by black dots 

Remark 1. U~ = Ui and = f/m+i are always admissible. In this case, the stationary states 
are the same as the corresponding initial states, and there are no waves. 

Remark 2. Out-flux qi = mm{DY ,S^} < D,- and in-flux qm+i = ^^^{^m+i^'^m+i} — "^m+i- 
hat is, D j is the maxi mum sending flow and S,n+i is the maximum receiving flow in the sense of 



(|Daganzo . 



1994 



1995h . 



Remark 3. The "invariance principle" in (ILebacque and Khoshyaran , 



20051) can be interpreted 



as follows: if D- = Q, then q{U^ ) < Dj; if Sj^^j = C^+i, then l{U^^^) < 5,-. We can see that 
Theorem 12. H is consistent with the "invariance principle". 



Corollary 2.2 For the upstream link i, qi < Di,- qi < Di if and only ifU- = {Ci.qi), and qi = Di 
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if and only ifU- = {Di.Ci). For the downstream link m+ \, qm+i < Sm+i; <?m+l < 5'ra+l if and 
only ifU^^^ = (<?„+!, C„+i), and q^+i = S,n+i if and only ifU^^^ = (Cm+i,5,„+i). That is, given 
out-fluxes and in-fluxes, the stationary states can be uniquely determined. 

For interior states, the waves of the Riemann problem on link i with left and right initial con- 
ditions of t/-^ and Ui{0^,t) cannot have negative speeds. Similarly, the waves of the Riemann 
problem on link m+l with left and right initial conditions of U,„+iiO^ and U^^^ cannot have 
positive speeds. Therefore, interior states t//(0^,?) and Um+i{0^,t) should satisfy the following 
admissible conditions. 

Theorem 2.3 (Admissible interior states) For stationary states Uf' and U^_^^, interior states 
Ui{0^ ,t) and (0+,?) are admissible if and only if 



The proof is quite straightforward and omitted here. The regions of admissible upstream interior 
states in both supply-demand and fundamental diagrams are shown in Figure [6l and the regions 
of admissible downstream interior states are shown in Figure Ul From the figures, we can also 
determine the types and traveling directions of waves with given stationary and interior states on 
all links, but these waves are suppressed and cannot be observed. Note that, however, we are able 
to observe possible interior states in numerical solutions. 

Remark 1. Note that Ui{Q^ ,t) = t/-^ and t/,„+i(0+,?) = U'^^^ are always admissible. In this 
case, the interior states are the same as the stationary states. 




(C,-,S- ) = Uj , when U- is SOC; i.e., S- < D- =C 

{Di (Q- , , 5,(0" , ) , when Ur is UC; i. e., Dr < 5r = Q 



(11) 



where Si{0 ,t)>D- , and 




where Dm_|_i(0+,?) > 



'm+l • 
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Si 

Ci 



nU,{0-,t) 




nU- = U,{0-,t) 



Ci D 




(C) 



(d) 



Figure 6: Admissible interior states for upstream link /: marked by black dots 

2.2 Entropy conditions consistent with discrete merge models 

In order to uniquely determine the solutions of stationary states, we introduce a so-called entropy 
condition in interior states as follows: 



q, = F,{Ui{0-j),---MO~.t).U,n+i{0+j)). 



(13) 



The entropy condition can also be written as 



q^+i = F{UiiO-,t),---MO-,t),U^+i{0+,t)), 

<? ; Ot;^ m+ 1 5 
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C„.+ i D 



(a) 



s 



(c) 



A„+i 



(b) 




(d) 



Figure 7: Admissible interior states for downstream link m + 1 : marked by black dots 

where a, e [0, 1] and J^^j a,- = 1 . We can see that the entropy condition uses "local" information in 
the sense that it determines boundary fluxes from interior states. In the discrete merge model, the 
entropy condition is used to determine boundary fluxes from cells contingent to the merging junc- 
tion. Thus, i^(i7i(0^,?),---,t/m(0^,f),i7„,+i(0+,?)) in (fT3l) can be considered as local, discrete 
flux functions. 



Here the boundary fluxes can be obtained with existing merge models. In (iDaganzol Il995h . 

F{U\{0^ ^t)^- ■ ■ ,Um{0^ ,t),Um+\{0^ ,t)) was proposed to solve the following local optimization 
problem 



max {qm+i} 

f^r.f^m+l'f^'-(0-,f),C/m+l(0+,0 



(14) 
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subject to 



Thus, we obtain the total flux as 

m 

F(f/i(0-,0,---,f/™(0",0,f/«.+i(0+,0) = min{£A-(0-,0,S™+i(0+,0}. 



;=1 



The proportions, «/, can be determ ined by so-ca 



total flux to each upstre am link. In dNewell . 



its m aximum flow. In (|Daganzd . 



led "distribution schemes", which distribute the 



19931), an on-ramp is given total priority an d can senc 



19951) . a priority -based scheme was proposed. In (|Lebacque . 



1996|) . a general scheme was proposed, a nd it was sugg ested to distribute out-fluxes according to 



the number of lanes of upstream links. In (IBanks 



2000f). an on-ramp is given total priority, but its 



20031), a fair sch eme is proposed 



flow is also restricted by the metering rate. In (|Jin and Zhangl . 
to dis tribute out-fluxes according to local traffic demands of upstream links. In (|Ni and Leonard . 
20051), a fair share of the downstream supply is assigned to each upstream proportional to its capac- 
ity, and out- fluxes are then determined by comparing the corresponding fair shares and demands. 
With these distribution schemes, the entropy condition, (fT3l) . yields unique solutions of boundary 
fluxes with given interior states. 



2.3 Summary of the solution framework 

To solve the Riemann problem with the initial conditions in ©-([T]), we will first find stationary and 
interior states that satisfy the aforementioned entropy condition, admissible conditions, and traffic 
conservation equations. Then the kinematic wave on each link will be determined by the Riemann 
problem of the corresponding LWR model with initial and stationary states as initial conditions 



( Jin et al. 



20091) . Here we will only focus on solving the stationary states on all links. We can see 



that the feasible domains of stationary and interior states are independent of the upstream supply. 
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Si, and the downstream demand, D,n+\. That is, the same upstream demand and downstream 
supply will yield the same solutions of stationary and interior states. However, the upstream and 
downstream wave types and speeds on each can be related to 5/ as shown in Figure l^d) and D^+i 
as shown in Figure |2d). 



3 Solutions of two merge models 

In this paper, we solve the Riemann problem for a merging junctions with two upstream links; i.e., 
m = 2. Different merge models have different entropy conditions (fT3] ). In this section, we consider 
two entropy conditions, i.e., two merge models. Here we attempt to find the relationships between 
the boundary fluxes and the initial conditions. 

qi=Fi{UuU2,U3). (15) 

In contrast to local, discrete flux functions i^(t/i(0^,?), 6^2(0 ^O? ^3(0^,0)' FiiUi,U2, C/3) can be 
considered as global, continuous. With the global, continuous fluxes, we can find stationary states 
from Corollary 12.21 With the solution framework in the preceding section, we can then find the 
kinematic waves of the Riemann problem with initial conditions {Ui.Uj, U3). 



3.1 A fair merge model 



We consider the fair merging rule proposed in (jjin and ZhangL |2003|) . in which 



i= 1,2. 



(16) 



z)i(o-,o+£>2(o-,0 

Obviously, the optimization problem of (fT3l) subject to the fair merging rule yields the following 



solutions 



qi = min{ 1 , 



53(0+,0 



z)i(o-,o+£>2(o-,0 



}A(0-,0, i=h2. 



(17) 
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Thus in the Riemann solutions, stationary and interior states have to satisfy (|T7l) . traffic conserva- 
tion, and the corresponding admissible conditions. 

Theorem 3.1 For the Riemann problem with the initial conditions in ® and ([7|), stationary and 
interior states satisfying the entropy condition in (123), traffic conservation equations, and the cor- 
responding admissible conditions are the following: 

1. When Di+D2< S3, Ur = Ui{0-j) = {Dud) (i = l,2j and = UsiO+j) = {Di + 
D2.C3); 

2. When Di +D2 = S3, Ur = U,{0-j) = {D,,Q) (i = \,2), U+ = (Cs^Ss), UsiO+j) = (Cs^Ss) 
or (1)3(0+,?), 53(0+,?)) withD^{Q+,t) > S3 andS^iQ+^t) > S3 when S3 < C3; 

3. WhenD, > ^^3 (/ =1,2), U+ = C/3(0+,0 = (C3,S3), andUr = U,{0~j) = (Q, cjtq'^s)- 

4. When Di + D2 > 53 and Di < ^^-^3 (^ j = I or 2 and i ^ j), = C/3(0+, t) = (C3, S3), 
ur = {D,,Q), f//(0-,0 = is^p^Q), and UJ = Uj{Q-,t) = {Cj,S3-Di). 

The proof of the theorem is given in Appendix A. In Figure [8l we demonstrate how to obtain the 
pair of (Di,D2) from initial states Ui and U2 and how to obtain the line of 53 — 53 in the supply- 
demand diagrams. According to the relationship between {D\,D2) and ^3, the domain of (D 1,^2) 
can be divided into four regions. In FigureHl we further demonstrate how to ob t ain sta tionary states 
for all these four regions. The results are consistent with those in (IDaganzoll 19961) in principle: 
for initial conditions in region I, there are no queues on upstream links; for initial states in region 
II, there is a queue on link 2; for initial states in region III, there are queues on both links; and, for 
initial states in region IV, there is a queue on link 1. 

Corollary 3.2 For the Riemann problem with initial conditions in ® and ([71), the boundary fluxes 
satisfying the entropy condition in 4771) . traffic conservation equations, and the corresponding 
admissible conditions are in the following: 
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Figure 8: Four types of solutions for a fair merge model 



1. When Di +D2 < S^, qi = Di (i = l,2j and = Di +D2; 

2. When Di > c^^^ = 1 , 2j, <?/ = cl^-^s and q^ = S3; 

3. When Di +D2 > ^3 and Di < (/, j = 1 or 2 and i 7^ j), qi = Di, qj = S3 — Di, and 
q3 = S3- 



That is, for /, j=\or2 and i 7^ j. 



min{D,, maxjS's — Dj, 



S3 



-Ci]]. 



(18) 



C1+C2 

The solutions of fluxes in four different regions are shown in Figure [lOl in which the dots represent 
the initial conditions in (Di,D2), and the end points of the arrows represent the solutions of fluxes 



{qi,q2). The fluxes in (fT8l) are exactly the same as in (INi and Leonard. 



2005h. In this sense, the 



discrete merge model with (flTl) converges to the merge model in (|Ni and Leonard . 
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[x+ = (C3.a,) 



(C) 



(d) 



Figure 9: Solutions of stationary states for different initial conditions for the fair merge model 



Corollary 3.3 If Uj (i = 1,2) and satisfy 



S3 



mm{Di,Si} = min{A,max{53-D^,— — — Q}}, 

Ci +L2 

S3 S3 
min{Z)3,S3} = min{Di,max{53 — D2, — — Ci}} + min{D2,max{53 — Di, ;;;; — l-;^^}} 



C1+C2 



C1+C2 



min{Di +^2,53}, 



then the unique stationary states are the same as the initial states, and traffic dynamics at the 
merging junction are stationary. 
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Figure 10: Solutions of fluxes for the fair merge model 



Proof. The results follow from Theorem 13. 1[ from which we can see that stationary states U- and 
5^ satisfy 



"^3 



min{D- ,5,- } = min{D- ,max{5;^ -D ■ , — ^— C,}} 
min{D^,5^} = min{Dj^,max{5-^ — 



5+ 5+ 
3 Ci}} + min{D2,max{5^-Df,^^^^C2}} 



C1+C2 



min{D7+D2 ,5^}. 



C1+C2 



3.2 A constant merge model 



We consider a merge model proposed in (ILebacqud . 



19961), in which 



^, = min{A(0",0,a,:S3(0+,0}, /= 1,2 
where a, are constant distribution proportions, a, G [0, 1], and a, = 1. 



(19) 
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Theorem 3.4 For the Riemann problem with the initial conditions in (0) and (0), boundary fluxes 
satisfying the entropy condition in 4791) . traffic conservation equations, and the corresponding 
admissible conditions are the following: 

1. When D1+D2 < S3 and Di < ttiCi, (i= \ ,2), qi = Di and = Di+ D2; 

2. When Di > ttiCj, and Dj < — ttiC^ (i,j = 1 or 2 and i 7^ j), qt = CXiCj,, qj = Dj, and 
qs = UjCs+Dj. 

3. When Di + D2 > S3, St, — CtjCs < Di < ttiSs ( z, j =1 or 2 and i 7^ j), qi = Di, qj = S3— Di, 
and q3 = S3. 

4. When Di > tt/Ss (i = 1, 2j, qi = (Z/Ss, and ^3 = 53. 
That is, for /, j = \ or 2 and i 7^ j, 

qi = min{A-,aA,max{53-D^-,a,-53}}. (20) 

The proof of the theorem is in Appendix B. The solutions of fluxes in six different regions are 
shown in Figure [TTl in which the starting point of an arrow represents initial conditions, {Di.Dj), 
and the ending point represents the solution of fluxes, {qi.qi). With the fluxes, we can easily 
find all stationary and interior states as in Theorem |3.1[ Compared with the fair merging rule, the 
constant merging rule yield sub-optimal solutions in regions II and VI in Figure [TTl in which q3 < 
min{Di +^2,53}. That is, this merge model does not satisfy the optimization entropy condition 
in ([111). 

4 Invariant merge models 

For the two merge models in the preceding section, the local, discrete flux functions ■, ■) are 
not the same as the global, continuous flux functions - , ■). That is, boundary fluxes computed 
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Figure 1 1 : Solutions of fluxes for the constant merge model 

from local flux functions are not constant. In this sense, the fair and constant merge models are not 
invariant. In the following, we propose two invariant merge models, in which ■) and Fi{-, ■, ■) 
have the same functional form. 



4.1 An optimal invariant merge model 



We consider the following priority-based merge model (.Daganzo . 
1 , or 2 and / ^ j 



1995 



19961), in which i,j 



^, = min{A(0",f),max{53(0+,0 -D^Q-,?), aA(0+,0}} 



(21) 



where a, are priority distribution proportions, a,- G [0, 1], and a, = 1. 



Theorem 4.1 For the Riemann problem with the initial conditions in (0) and (l7|), boundary fluxes 
satisfying the entropy condition in f [27]) . the traffic conservation equations, and the corresponding 
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admissible conditions are the following: 



qi = min{D/,max{53 —Dj, tt/Ss}}, ij = 1, or 2, and i ^ j 



(22) 



The proof of the theorem is given in Appendix C. Since - , ■) = - , ■), the merge model 
(|2TI) is invariant. From the proof, we can see that the merging rule is optimal in the sense that 
^3 = min{Di +^2,53}. The solutions of fluxes in four different regions are shown in Figure 
[T2I in which the starting point of an arrow represents the initial demands, {Di.Di), and the end 
point represents the solutions of fluxes, {qi,q2). Note that, when a, = r,'^'c2 ~ ^ ^ .^^^^ . ^^^ 



distribution scheme is the same as (fT8l) . Therefore, the merge model in (iNi and Leonard . 
invariant. 



2005) is 




C?, Di,qi 



Figure 12: Solutions of fluxes for the priority -based merge model 



When a, > q+c^' merge model (|2TI) gives higher priority to vehicles fro m upstream link 
Thus, (I2T]) is another representation of the priority merg e model 



19961). For an extreme merge model with tti = 1 and 0^2 = (.Newelli 



jroposed in ( Daganzd . 



1993 



BanksL 



1995L 



20001), vehicles 
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from link 1 is given absolute priorities. In this case, we have 

<7i = min{Di,53}, 
q2 = min{Z)2,max{53 -Di,0}}. 

4.2 A suboptimal invariant merge model 

Similarly, the continuous version of the constant merge model will lead to another invariant merge 
model, in which 

^, = min{A(0",0,aA,max{53(0+,0-^y(0",0,«A(0+,0}}, = 1, or 2, mdi^j(23) 
with tti e [0, 1], and = 1. 

Theorem 4.2 For the Riemann problem with initial conditions in (0) and (0), the boundary fluxes 
satisfying the entropy condition in ([27]), traffic conservation equations, and the corresponding 
admissible conditions are in the following: 

qi = min{D,-, a,C3, maxjSs — Dj, tt/Ss}}, j = 1 , or 2, and i ^ j (24) 

The proof of the theorem is omitted. The merge model (l23l) is obviously invariant. In addition, 

the solutions of fluxes in six different regions are the same as in Figure \TT\ We can see that this 

merge model is suboptimal. Assuming that both links 1 and 3 are the mainline freeway with 

capacity Ci = C3 and Di =83 = Ci, and link 2 is a metered on-ramp with a metering rate D2 < 

C2 < Ci. When Di > OCiCs, and the metering rate D2 < 0C2Cti, then qi = aiC^ and q2 = D2. In this 

case, q3 = CtiCs +D2, and the utilization rate of link 3's capacity is 

qs _ «iC3+Z)2 _ aiC3+D2 ^ 

min{Di+D2,C3} ~ min{Di +D2,C3} ~ C3 " 

That is, as much as a2 of link 3's capacity is wasted. In contrast, if we increase the metering 
rate of link 2, such that D2 > a2Cj,. Then (l23l) predicts ^3 = C3. That is, low metering rates 
could cause a lower utilization rate of link 3's capacity. This property of the merge model (|23l) is 
counter-intuitive, and the merge model is not physical. 
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5 Numerical examples 

In this section, we numerically solve various merge models and demonstrate the validity of our 
analytical results. Here, both links 1 and 3 are two-lane mainline freeways with a corres ponding 



normalized maximum sensitivity fundamental diagram (|Del Castillo and Benitea . 

[0,2]) 



1995h is (p e 



Qip) 



p < 1 — exp 



Link 2 is a one-lane on-ramp with a fundamental diagram as (p G [0, 1]) 

1 



QiP) 



p < 1 — exp 



,1,1 
l-exp(-(--l. 



Note that here the free flow speed on the on-ramp is half of that on the mainline freeway, which is 
1. Thus we have the capacities Ci = C3 = 4C2 = 0.3365 and the corresponding critical densities 
pel = Pc3 = 2pc2 = 0.4876. The length of all three links is the same as L = 10, and the simulation 
time is r = 360. 

In the following numerical examples, we discretize each link into M cells and divide the sim- 
ulation time interval T into steps. The time step At — T /N and the cell size Ax = L/M, with 



At = 0.9 Ax, satisfy the CFL condition (|Courant et al. 



1928b 



At At 
Ax Ax 



0.9 < 1. 



Then we use the following finite difference equation for link / = 1,2,3: 



Pi,m + ^ 1 /2 ~ ^1,(11+1/21 



where p"^^^ is the average density in cell m of link / at time step n, and the boundary fluxes ^"„j_i/2 
are determined by supply-demand methods. For example, for upstream links i =1 and 2, the in- 
fluxes are 
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where is the demand of cell m — 1 on link /, 5"^ the supply of cell m, and D"q the demand 

of link /. For link 3, the out- fluxes are 

?3,m+i/2 = min{D^^„,5^^„,+i}, m=l,---,M, 

where S^j^^^ is the supply of the destination. Then the out-fluxes of the upstream links and the 
in-flux of the downstream link are determined by merge models, which are discrete versions of 

^'i,M+l/2 = '^■(^l,M'-^2,M''^3j)' 
^3,1/2 = '?l,M+l/2 +^2,M+l/2- 

In our numerical studies, we only consider the fair merge model (fTTI) and its invariant counter- 
part (|2T)) . That is, in the fair merge model, we have 

D" 

<M+i/2 = min{Z);^^ „ 

^3,1/2 = ^^{D'Im+DIm.SIi}. 

In the invariant fair merge model, we have (/, j = 1,2 and / ^ j) 

<M+i/2 = min{D;'^,max{5^ i-D;;^,a,5^ J}, 
^3,1/2 = ^HD'Im+DIm,SIi}, 

where a, = C// (Ci + C2) . 

5.1 Kinematic waves, stationary states, and interior states in the fair merge 
model 

In this subsection, we study numerical solutions of the fair merge model in (flTl) . Initially, links 1 
and 3 carry UC flows with pi = P3 = 0.35, and at ? = a traffic stream on link 2 starts to merge 
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into link 3 with p2 = 0.1. That is, the initial conditions in supply-demand space is Ui = Ut, 
(0.3131,0.3365) and t/?. = (0.0500,0.0841'). H ere we use the Neumann boundary condition in 



supply and demand (|Colella and Pucketi 



9004 V — D" D" — n" and <>" — <>" 



Therefore, we have a Riemann problem here. 

In this case, Di +D2 = 0.3631 > ^3 = 0.3365 and < c^^3 = 0-0673. Thus according 
to Theorem 13. 1[ we have the following stationary and interior states = Uj,{0^,t) = {€3,83) = 
(Ci,Ci), = (D2,C2), t/2(0-,0 = (5^1)2, C2), and t/f = t/i(0-,0 = {CUS3-D2). Further- 
more, from the LWR model, there will be a back-traveling shock wave on link 1 connecting Ui to 
i7f , no wave on link 2, and a rarefaction wave on link 3 connecting to U3. 



(a) Contour plot of p 




(c) Contour plot of p 



(b) Contour plot of p- 





Figure 13: Solutions of the fair merge model ([IT]): M = 160, = 6400. 



In Figure [T3l the solutions of pi, p2, and p3 are demonstrated with M = 160 and = 6400. 
From the figures, we can clearly observe the predicted kinematic waves. In Figure \T3th) there is 
a very thin layer of higher densities in the last cell of the upstream link near the merging junc- 



27 



tion, which is caused by the interior state as predicted. In addition, we can observe at t = T 
the approximate asymptotic values: U[ = Ui{0^ ,t) = (0.3365,0.2865), and = pi(0^,?) — 
0.8277; = (0.05,0.0841), U2{0-,t) = (0.0587,0.0841), = 0.1, p2(0",0 = 0.1179; U+ = 
(0.3365,0.3365), and = pc3 = 0.4874. These numbers are all very close to the theoretical val- 
ues and get closer if we reduce Ax or increase T. That is, the results are consistent with theoretical 
results asymptotically. 



a 0.048 7 



IT 



(a) Out-flux of link 2 



IT*' 



Ax=1 

Ax=1/4 

Ax=1/16 



10 
t 



(b) Density in the downstream cell of link 





Ax=1 

Ax=1/4 

Ax=1/16 



10 
t 



Figure 14: Evolution of the out-flux and the density in the downstream cell of link 2 for the fair 
merge model (fTT] ) 



In Figure [TH we demonstrate the evolution of the out-flux and the density in the downstream 
cell of link 2 for three different cell sizes. From Figure H^la) we can see that, initially, the out-flux 
of link 2 is 

-53 = 0.3365 = 0.0463, 



D1+D2 0.05 + 0.3131 
which is not the same but approaches the asymptotic out-flux D2 = 0.05. Correspondingly the 
density in the downstream cell of link 2 approaches the interior state, as shown in Figure [T4T b). 
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Furthermore, as we decrease the cell size, the numerical results are closer to the theoretical ones at 
the same time instant. This figure shows that the fair merge model is not invariant, but approaches 
its invariant counterpart asymptotically. Note that the densities in any other cells of link 2 remain 
constant at 0.1. 

5.2 Comparison of non-invariant and invariant merge models 

In this subsection, we compare the numerical solutions of the fair merge model (flTI) with its in- 
variant counterpart ((22)) . Initially, links 1 and 3 carry UC flows with pi = Ps = 0.35, and at 
t = a traffic stream on link 2 starts to merge into link 3 with p2 = 0.1. Different from the ex- 
ample in the preceding subsection, here we use the following boundary conditions: D" q = ^, 
^2 = 0.05 -I- 0.03 sin(n;rA?/60), and ^^j^-^ = p^. Thus we have a periodic demand on link 2. 

0.12 
0.1 
0.08 
^ 0.06 
0.04 
0.02 


50 100 150 200 250 300 350 

t 

Figure 15: Difference in the solutions between the fair merge model (fTTI) and its invariant counter- 
part (EU) 

We use pf„j for the discrete density from the fair merge model (flTl) and p"^ from its invariant 
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counterpart (|22l) . Then we denote the difference between the two solutions by 

3 M 

£{nM) = Y.ll\Pln-Pun\^- (25) 

i=l m=\ 

In Figure [151 we can see that the difference decreases if we decreases the cell size. This clearly 
demonstrates that the fair merge model (fTTI) converges to its invariant counterpart (|2TI) . 

6 Conclusion 

In this paper, we studied continuous kinematic wave models of merging traffic flow which are 
consistent with discrete CTM with various distribution schemes. In particular, we introduced the 
supply-demand diagram of traffic flow and proposed a solution framework for the Riemann prob- 
lem of merging traffic flow. In the Riemann solutions, each link can have two new states, an interior 
state and a stationary state, and the kinematic waves on a link are determined by the initial state and 
the stationary state. We then derived admissible conditions for interior and stationary states and 
introduced various distribution schemes as entropy conditions defined in the interior states. Then 
we proved that stationary states and boundary fluxes exist and are unique for the Riemann problem 
for both the fair and constant distribution schemes. We also discussed two invariant merge models, 
in which the local and discrete flux is the same as the global and continuous flux. With numeri- 
cal examples, we demonstrated the validity of the proposed analytical framework and that the fair 
merge model converges to its invariant counterpart when we decrease the cell size. Compared with 
existing discrete kinematic wave models (i.e. CTM with various distribution schemes) of merging 
traffic flow, the continuous models can provide analytical insights on kinematic waves arising at 
merging junctions; and, compared with existing analytical models, they are physically meaningful 
and consistent with existing CTM. 

In this study, we introduced a new definition of invariant merge models, in which fluxes com- 
puted by discrete models should be the same as those by their continuous counterparts. An impor- 
tant observation is that both the fair and constant merge models are not invariant. For example, for 
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the fair merge model, at ? = the local fluxes from (flTl) are 



qi = min{ 1 , 



^3 



}Di 



/= 1,2, 



which are different from (1181) when only one upstream is congested; i.e., when Di +D2 > S3 and 
Di < Ci+C2 '^3" However, the results here suggest that the discrete fluxes converge to the continuous 
ones after a sufficient amount of time or at a given time but with decreasing time intervals. That 
is, the non-invariant discrete merge models do not provide incorrect solutions but just approximate 
solutions to the corresponding continuous merge models. Compared with invariant discrete merge 
models, the fair merging model has some unique mer its; e.g., it can be easi l y exte nded to general 



junctions with multiple upstream links as shown in (|Jin and Zhan§ 



demonstrated in (|jin et al.L 



2001 



2004). Note that, as 



2009h . invariant models can also yield interior states. 
This paper presents a systematic framework for solving kinematic waves arising from merging 
traffic in supply-demand space. We expect that boundary fluxes, stationary states, and kinematic 
waves for other distribution schemes can also be solved in this framework. For example, we can 
obtain the stationary states and kinematic waves for the following artificial merge model 

0.953(0+,0 



qi = min{l, 



Di(0-,0+£>2(0-,? 

in which only 90% of the downstream supply can be utilized by the upstream traffic. 

The Riemann problem for a merge with three or more upstream links can be discussed in 
the same framework, but, due to the space limitations, this will be discussed in future studies. 
Generally, there are systematic lane-changing activities at merging junctions, and ex isting merge 



mode ls based on supply-demand method cannot capture the impacts of lane-changes (|Laval et al. 



20051). In the future, it would be interesting to analyze the formation and dissipation of traffic 



queues at merging junctions when lane-changes are considered. 
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Appendix A: Proof of Theorem I3J^ 



Proof. From traffic conservation equations in ([8]) and admissible conditions of stationary states, we 
can see that 

^3 <min{Z)i -1-^2,53}. 

We first demonstrate that it is not possible that q2, < min{Di + , 5*3 } < min{Ci + C2 , C3 } . Other- 
wise, from ([Jo]) and we have t/3(0+,?) = U:^ = (qs^C^) with < S3; Since q{U^)+q{U2) = 
q2, < D\ +D2, then we have q{U^) < D,- for at least one upstream link, e.g., q\ < D\. From ^ 
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and (fTTI) we have t/i (0 ,t) = = (Ci ,qi). Then from the entropy condition in (fTTI) we have 

^3 = min{Ci +D2(0",0,C3}, 
Q] — minjl, r}C\. 

Since qj, < C3, from the first equation we have ^3 = Ci +^2(0^,?) < C3, and from the second 
equation we have qi = Ci, which contradicts qi < Di. Therefore, 

<?i +<?2 = <?3 = min{Z)i +^2,53}. 

That is, the fair distribution scheme yields the optimal fluxes for any initial conditions. 

(1) When Di +D2 < S3, we have ^3 = £>i +£>2 < S3. From ^ and ^ we have t/3(0+,0 = 

= {Di +D2,C3). Since q\+q2= Di +D2 and qt < Di, we have qi = Di. From ^ we 
have Ur = (A-,C,). From (EJ) we have Ui{Q- ,t) = {Di{0- j),Si{0- ,t)) with Si{0-,t) > 
DJ = Di. From (fTTI) we have 

q3 = mm{Di{0-,t)+D2{0-,t),C3}=Di+D2<S3<C3, 

* = "'""■ D.(0-.,)+D.(0-,» '°-'°''^'="'- 

Thus, A(0",0 = A- < Si{0',t). Then t//(0",0 = Ur = {Di.Q). In this case, there are no 
interior states on all links. 

(2) When Di +D2 = ^3, we have ^3 = ^3. From ^ and CH), we have = (C3,53) and 
U3{0+j) = (D3(0+,0,53(0+,0) withD3(0+,0 >5+ =53. Since ^1 +^2 = ^1 +^2 and 
qi < Di, we have qi = Di. From Q and (fTTI) . we have Uf' = (Di^Q) and Ui{Q^ ,t) — 
{Di{Q-,t),Si{Q-,t)) with 5,(0-,f) > Dr = Di. From flT]) we have 

^3 = min{D^{0-,t)+D2iO-,t),S3iO+,t)} = D^+D2 = S3, 

* - D.(o-':r+4-.>) '^''°''''^°'- 

We can have the following two scenarios. 
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(i) li Di{0-,t)+D2{0-,t) > 53(0+,?) = Di +D2 = 53 < D3{0+,t), then U3{0+,t) = 

= (C3, 53) and there is no interior state on link 3. Moreover, we have 

whichleadstoA(0~,?) < A<5',(0",?). From the assumption that Di (0", +£>2 (0^,0 > 
Di +D2, we have A(0",0 = Z),-. Further we have t/,(0",0 = Uf' = {DuQ), and there 
are no interior states on links 1 or 2. 

(ii) If 53(0+, > Di(0-,f) +D2(0-,0 = Di +D2 = 53, A(0-,0 = Thus f/,(0^,0 = 
Uj^ = (Di.Ci), and there are no interior states on links 1 or 2. Moreover, t/3(0+,?) 
satisfies 53(0+,?) > 53 and ^3(0+,?) > 53. Thus there can be multiple interior states 
on link 3 when 53 < C3. 

(3,4) When D[+D2> S3, for upstream links, at least one of the stationary states is SOC. Other- 
wise, from dll) we have t/-^ = (A,Q), and qi+q2 = Di + A > S3, which is impossible. In 
addition, we have ^3 = 53 < Di +^2. From (flOl) we have = (€3,83). From (fT2l) we have 
[/3(0+,f) = (A(0+,0,'53(0+,0) with A(0+,0 >'5^ = '53- From (di]) we have 

q3 = mm{Di{0-,t)+D2{0-,t),S3{0+j)} = S3<Di+D2, 

- D.(0-T+4-.>) '^''°''''- 

If ^1(0-,?)+ A(0",0 < 53(0+,?), then Di(0-,f)+ A(0",0 = ^3 < Di +A and = 
A(0^,0- This is not possible for the SOC stationary state Uf' = Uj{0^,t) = {Ci,qi) with^; < 
A < Q. Thus 53(0+,?) < Di(0-,0 + A(0",0. S3{0+,t) = S3<Di + A, and U3{0+j) = 
i/j'" = (C3, 53) . Hence for both upstream links 

«■ = D.(0-,,) + D,(0-,0'''<°"'''- 
(3) When D, > q^-^^S3 (i = 1,2), stationary states on both links 1 and 2 are SOC with 
U~ = Ui{0^,t) = (Ci.qi) with qi < Di. Otherwise, we assume that link 1 is SOC with 
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Ui{0-,t) = C/f = (Ci,^i) and link 2 is UC with = (^2,C2). Then 

Di = 7^2(0 ,0 < — C2 < D2, 

Ci+D2(0-,0 ' ^ - C1+C2 

which is impossible. From (fTTI) . we have 

^'"Ci+C2^" 

andC//(O-,O = f/r = (Q,^0- 
(4) When Di + D2 > S3 and D,- < q^^^S^ = 1 or 2 and / 7^ j), we can show that 
stationary states on links j and / are SOC and UC respectively with Uj' = Uj{0^j) = 
{Cj,qj) with qj < Dj, Ur = [DuQ], and 5,(0^,0 > A'. Otherwise, t/,(0",0 = Ur = 
{Ci,qi) with qi < Di, and 

:Ci > 77^53 > A, 



G + C1+C2 

which is impossible. Since at least one of the upstream links has SOC stationary state, 
the stationary states on links i and j are UC and SOC respectively. From (fTTI) . we have 
a unique interior state on link z, t/,(0^,?) = (^p^Cy,C/), and qj = S3 —Di. 

For the four cases, it is straightforward to show that (fTSi ) always holds. ■ 



Appendix B: Proof of Theorem 1X4 



Proof. 



(1) When Di +D2 < S3 and A < aA (z = 1,2), then q^ < S3, and U+ = f/3(0+,0 = (^3,C3). 
Assuming that < A < Q» then Uf' = A(0^,0 = (QiQi)- From (fT9l ). we have = 
min{ C,-, OCiCs} = a^Cs < Di, which contradicts Di < a,C3. Therefore, qi = Di and q3 = 
D1+D2. 
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(2) When D/ > a,C3 and Dj < 53 — a^Cs (/, 7 = 1 or 2 and / 7^ 7), we first prove that qi < Di and 
then qj = Dj. 

(i) If qi = Di, then from we have tt/Cs < D/ = min{D/ (0" , ? ) , tt/Ss (0+ , ? ) } < ttiSs (0+ , ? ) 
which leads to C3 < 53(0+,?). This contradicts C3 > 53(0+,?). Thus, ^/ < D,, and 
f/r = [/,(0-,0 = (Q,^,). 

(ii) If qj < Dj < Cj, then UJ' = Uj (0" , r ) = {Cj, qj). From ([l9l), we have qj = min{Q, ajS-i (0 
a^-53(0+,0 < D^- < 53 - a,C3 < ayC3. Thus 53(0+, < C3. Since = a,-53(0+,0, we 
have <5r3 = qi + qj < a,:53(0+,r) + ^3 - < S3, which leads to U^^ = [73(0+,?) = 
(^3,C3). This contradicts 53(0+,?) < C3. Therefore, qj =Dj. 

From ([HI), we have <7,- = a,-53(0+,f). Thus ^3 = + < ^3 - ocA + a/53(0+,0 < 53. 
Then, U+ = f/3(0+,0 = (^3,^3), qi = OcA, qs = CCiCs +Dj. 

(3) When Z)i +D2 > ^3, S3 - ajC3 < Di < aiSs (ij =1 or 2 and i ^ j), we first prove that 
q^ = Si and then that qi = Di. Therefore qi = Di, and qj = S^— Di. 

(i) If qT, < 53, then = U2,{0^j) = {q^.CT,). We first prove that at least one upstream 
stationary state is SOC and then that none of the upstream stationary states can be SOC. 
Therefore, ^3 = 53. 

(a) If none of the upstream stationary states are SOC, then qi = Di (z = 1,2), and 
<13 = Qi + Q2 = Di + D2 > St,, which contradicts ^3 < 53. 

(b) Assuming that qi < Di, then Uf' = t/,(0^,?) = {Ci.qi). From (fT9l) . we have qi = 
min{C,, OCiC-i] = UiCj, < Di. This is not possible, since D, < UiSj,. Thus, qi = Di. 

(c) Assuming that qj < Dj, then Uj' = Uj(0^ j) = (Cj.qj). From (fT9l) . we have qj = 
min{Cj, CCjC^} = ajC^. Then Di = qi = ^3 — CXjC^ < 'S'3 — CXjC^, which contradicts 
A>53-a,C3. 
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(ii) Now assume^/ < A, then = Ui{0^,t) = {Q.qi). From (fT9l). we have q; = min{C,-. a,5^ (0+ . t) } 
oc,53(0+,?) < Di < UiS^. Therefore, qj < aj53(0+,?) < ajS^, which leads to qi + qj < 
53. This is not possible, since qi + qj = St,. 

(4) When D, > (XiS^, we first prove that q^ = and then that qi > a, 53 (z = 1,2). Therefore, 
qi = aiSs (z = 1 , 2), and <?3 = 'S'3. 

(i) If q2 < 53, then = U2,{0^,t) = (q^.C^). We first prove that at least one upstream 
stationary state is SOC and then that none of the upstream stationary states can be SOC. 
Therefore qi+q2 = q3 = S3. 

(a) If none of the upstream stationary states are SOC, then qi = Di (z = 1,2), and 
q3= Di+D2> S3, which contradicts ^3 < 53. 

(b) Assuming that qi < Di, then Uf' = Ui{0^,t) = {Q.qi). From (fT9l) . we have qi = 
min{Q, oc,C3} = a,C3. Then qj = ^3 — qi < S3 — a,C3 < ajS3 < Dj. Thus, UJ' = 
Uj{Q^ ,t) = {Cj,qj), and from (fT9l) we have qj = min{Cj, OCJC3} = ajC3 > ajS3, 
which contradicts qj < ajS3. Thus, t/^ is UC and qi = Di. 

(c) Similarly we can prove that Uj' is UC and qj = Dj. 

(ii) If qi < a/53 < Di < Q, then Ur = Ui{0' j) = (Q^qi). From we have qi = 
min{C/, a,53(0+,?)} = a/53(0+,?) < a/53. Therefore, qj < a^-53(0+,?) < ajS3, which 
leads to qi + qj < S3. This is not possible, since qi + qj = S3. 



For the four cases, it is straightforward to show that (|20l) always holds 



Appendix C: Proof of Theorem 14^1 

Proof. 
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First (I2TI) implies that 

^3 = mm{Di{0-,t)+D2{0-,t),S3{0+,t)}, 

which can be shown for three cases: (i)Di(0",f)+D2(0",?) < 53(0+,?), (ii)A(O",?) > a,-53(0+,?), 
and(iii)Z)i(0-,?)+/)2(0",0 >'^3(0+,0 andA(0",?) < aiS3{0+,t). 

(1) When Di +D2 < ^3, = qi + qi < Di + D2 < S3 < C3. Thus the downstream stationary 
state is SUC with = 11^ (0+ , f ) = (^3 , C3 ) . In the following, we prove that qi = Di, which 
is consistent with (|22l) . Therefore (|22|) is correct in this case. 

(i) Assuming that qi < Dj < Cj, then the stationary state on link i is SOC with = 
Ui{0^,t) = {Ci,qi). From (|2T1) . we have 

qi = min{C/, max{C3 — Dj (0^ , ? ) , UiC^ } } = max{C3 — Dj (0^ , t) , tt/Cs } < C,-, 
= min{Dy(0",0,max{C3-Q,ayC3}}. 

We show that the two equations have no solutions for either UjC^ < Dj (0^ ,t) or UjC^ > 
Dj{0-j). Thus qi = Di. 

(a) When ajC^ < Dj{0^,t), we have a,C3 > C3 —Dj{0^,t). From the first equation 
we have qt = UiC^. From the second equation we have qj = Dj{0^,t) > UjC^ 
or qj = max{C3 — C,-, CCjC^} > UjC^. Thus qi + qj > C3 > ^3, which contradicts 
q?, < S3. 

(b) When ajC3 > Dj(0^ j), we have a,C3 < C3 — D j{0^ j). From the first equation 
we have t// = C3 — Dy(0^,?). From the second equation we have qj = Dj{0^j). 
Thus qi + qj = C3, which contradicts ^3 < 53. 

(2) When D,- > a,-53, S3 — Dj < a/53. In the following we show that ^3 = 53 and = a, 53, 
which is consistent with (|22)) . 
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(i) If ^3 < S3, then the stationary state on link 3 is SUC with = Uj{0'^,t) = (q^^Cj). 
Also at least one of the upstream stationary states is SOC, since, otherwise, qi+q2 — 
Di+D2> 53. Here we assume Uf' = Ui{0^ j) = (Ci.qi). From (|2T)) we have 

qi = min{Q, max{C3 - Dj (0" , f ) , } } = max{C3 - Dj (0" , , aA } , 
qj = min{Z)^(0",f),max{C3 -C/,a^C3}}. 

We show that the two equations have no solutions for either C3 — Dj{0^j) > ttiCs or 
C3 -Dj{0-,t) < aA. Thus ^3 = 53. 

(a) If C3 —Dj{0^,t) > ttiC^, Dj{0^,t) < UjC^. From the first equation we have qt = 
C3 — Dj (0^ j). From the second equation we have qj = Dj (0^ , t) . Thus qi + qj = 
C3, which contradicts ^3 < ^3 < C3. 

(b) If C3 —Dj{Q^,t) < a,C3, Dj{Q^ ,t) > ajC^,. From the first equation we have qj = 
ttiC^ . From the second equation we have qj = Dj (0^ , > CU/Q or qj = max{C3 — 
Q, OCjC^} > ttjC^. Thus qi + qj > C3, which contradicts <?3 < ^3 < C3. 

(ii) If qi < a/53 < Di < C, for any « = 1,2, then Ur = [/,(0",?) = (Q,^,). From dH]) we 
have 

qi = max{53(0+,0-D/0-,0,a/53(0+,0}<Q, 
qj = min{DX0-,0,max{53(0+,0-Q,a;53(0+,f)}}. 

The first equation implies that a/53(0+,?) < ttiS^; i.e., 53(0+,?) < 53. In addition, 
53(0+, f) -D^(0",0 < a,-53. Thus, 53(0+,?) -C/ < ^3 -C,- < ^3 - ttiSs = ajS3, and 
max{5'3 (0+ , ? ) — C/, aj53 (0+ , } < 05^53 . From the second equation we have qj < ajSs . 
Thus qi + qj < S3, which contradicts qi + qj = S3. Thus qi > aiS3 for / = 1,2. Since 
qi + qj = S3, qi = a,-53. 

(3) When Di+Dj > S3 and D, < a, 53 for i,j= 1 or 2 and i ^ j. In the following we show that 
q3 = S3 and qi = Di, which is consistent with (|22|) . 
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(i) If ^3 < S3, then the stationary state on link 3 is SUC with = Uj{0'^,t) = {q^^Cj). 
We first prove that at least one upstream stationary state is SOC and then that none of 
the upstream stationary states can be SOC. Therefore, ^3 = 53. 

(a) If none of the upstream stationary states are SOC, then qi + q2 = Di + D2 > S^, 
which contradicts q^ < Sj. Thus, at least one of the upstream stationary states is 
SOC. 

(b) Assuming that qt < Di, then U^^ = Ui{Q^ ,t) = (Q^qi). From (|2TI) we have 

qi = min{C/,max{C3 -D;(0",?), OfA}} = max{C3 -Z)y(O",0, ocA} < A, 

which is not possible, since < ttiSi,. Thus qi = Di. 

(c) Assuming that (^r^- < D^-, then t/^^ = Uj{Q^ ,t) = (Cj^qj). Since^3 = min{D,(0^,?) + 
Cj.Ci} < St,, we have ^3 = D,(0^,?) +Cj < S3. From (|2TI) we have D,- = qi < 
Di{0-j). Thus Di + Dj < A(0",0 +Cj < S3, which contradicts Di+Dj > S3. 

(ii) If qi < Di, then Ur = t/,(0",?) = (Q,^,). From dlD, we have 

qi = max{53(0+,0 -£>;(0-,0,a/53(0+,0} < A < 05/53, 
qj = min{DX0-,0,max{53(0+,0-Q,a^-53(0+,0}}. 

From the first equation we have that 53(0+,?) < ^3. We show that the two equa- 
tions have no solutions for either 53(0+,?) — Dj{0^,t) > a/53(0+,?) or 53(0+, — 
Dj{0-,t) < a/53(0+,0- Therefore qi = Di. 

(a) WhenS3{0+,t)-Dj{0-,t) > a/53(0+,/'), we have D^- (0^,0 < ajS3{0+,t). Thus 
^- = 53(0+,?) -Dj{0-,t) and qj = Dj{0-,t). Then = 53(0+, < ^3, which 
contradicts qi + qj = S3. 

(b) When53(0+,0 -D/(0",?) < a/53(0+,0, wehave^,- = a/53(0+,f) and53(0+,0 - 
Ci < S3 (0+ , - qi = 0CjS3 (0+, t) . Thus qj < ajS3 (0+ , , and qi + qj < ^3 (0+ , < 
53, which contradicts ^3 =53. 
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